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1. The first problem relating to the motion of a solid body in a viscous liquid which 
was successfully attacked was that of a sphere, the solution of which was given by 
Professor Stokes in 1850, in his memoir ^^On the Effect of the Internal Friction of 
Fluids on Pendulums," ' Cambridge Phil. Soc. Trans.,' vol, 9, in the following cases : 
(i.) when the sphere is performing small oscillations along a straight line ; (ii.) when 
the sphere is constrained to move with uniform velocity in a straight line ; (iii.) 
when the sphere is surrounded by an infinite liquid and constrained to rotate with 
uniform angular velocity about a fixed diameter : it being supposed, in the last two 
cases, that sufficient time has elapsed for the motion to have become steady. In the 
same memoir he also discusses the motion of a cylinder and a disc. The same class 
of problems has also been considered by Meyer'"' and OBERBECK,t the latter of whom 
has obtained the solution in the case of the steady motion of an ellipsoid, which- 
moves parallel to any one of its principal axes with uniform velocity. The torsional 
oscillations about a fixed diameter, of a sphere which is either filled with liquid or is 
surrounded by an infinite liquid when slipping takes place at the surface of the sphere, 
forms the subject of a joint memoir by Helmholtz and Piotrowski.J 

Very little appears to have been effected with regard to the solution of problems 
in which a viscous liquid is set in motion in any given manner and then left to itself. 
The solution, when the liquid is bounded by a plane which moves parallel to itself, is 
given by Professor Stokes at the end of his memoir referred to above ; and the solu- 
tions of certain problems of two-dimensional motion have been given by Stea.rn.§ 
In the present paper I propose to obtain the solution for a sphere moving in a viscous 
liquid in the following cases : — (i.) when the sphere is moving in a straight line under 
the action of a constant force, such as gravity ; (ii.) when the sphere is surrounded by 
viscous liquid and is set in rotation about a fixed diameter and then left to itself. 

* ' Crelle, Journ. Math.,' vol. 73, p. 3L 
t ' Crelle, Joxirn. Math.,' vol. 81, p. 62. 
J ' Wissenscliaftl. Abliandl.,' vol. 1, p. 172. 
§ * Quart. Journ. Math.,' vol. 17, p. 90. 
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Throughout the present investigation terms involving the squares and products of 
the velocity will be neglected. This is of course not strictly justifiable, unless the 
velocity of the sphere is slov7 throughout the motion. If, therefore, the velocity is 
not slov7 the results obtained can only be regarded as a first approximation ; and a 
second approximation might be obtained by substituting the values of the component 
velocities hereafter obtained in the terms of the second order, and endeavouring to 
integrate the resulting equations. I do not, however, propose to consider this point 
in detail. 

2. In the first place it will be convenient to show that the equations of impulsive 
motion of a viscous liquid are the same as those of a perfect liquid. 

The general equations of motion of a viscous liquid are 

> ckt , dtt , die , dii ^y, 1 dp a ^ 

. dt dm dy dz p ax 

with two similar equations, where v is the kinematic coefiicient of viscosity. 

If we regard an impulsive force as the limit of a very large finite force which acts 
for a very short time r, and if we integrate the above equation between the limits 
T and 0, all the integrals will vanish except those in which the quantity to be inte- 
grated becomes infinite when r vanishes ; we thus obtain 

1 d l'^ 

u — Uc^A rr- pdr = 0. 

Putting f pdr = m where m is the impulsive pressure at any point of the liquid, 
we obtain 

H'TK 

p{u — %) + T~ = 0^ ^^*^ ^^'5 

which are the same equations as those which determine the impulsive pressure at any 

point of a perfect Hquid. 

3. Let us now suppose that a sphere of radius a, is surrounded by a viscous Hquid 
which is initially at rest, and let the sphere be constrained to move with uniform 
velocity V, in a straight hne. If the squares and products of the velocity of the 
liquid are neglected, Professor Stokes has shown that the current function rjj must 
satisfy the differential equation 



vi^here 



d^ ^ mi6 d / ^ d\ 

dr^ r^ d0\ dd 



and (r, 0) are polar coordinates of a point referred to the centre of the sphere as 
oaigin. 
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Let R, @ be the component velocities of the Hquid along and perpendicular to the 
radius vector ; then, if we assume that no slipping takes place at the surface of the 
sphere, the surface conditions are 

j^ — ~~-~-^^- --^ = V cos ^, (2) 

1 d'ylr 

@ -- _ — ^ — -7-- == — • V sin ^ (3) 

a sm dr ^ ^ 

Also, at infinity R and % must both vanish. 
These equations can be satisfied by putting 

\^ = (\^j -[- ^^ sin^ d^ (4) 

where i/^j and \}i^ are functions of r and t, which respectively satisfy the equations 

79 9 "~~" ^j • • • • • .. . . . \ D ) 

dr^ T^ V at ^ ^ 

The proper solution of (5) is t/*i =f(t)/r, which it will be convenient to write in the 
form 

^^ = ^^7"^^ J ^ X (^) ^XP- ( - ^V^^t) da, ..... (7) 

where x {^) is an arbitrary function, which will hereafter be determined. 

In order to obtain the solution of (6), let us put i/r^ = re'^^^^^'^dw/dr, where iv is a 
function of r alone ; substituting in (6), and integrating, we obtain 

rw = A cos X (r — a + ^)j 

where a is the radius of the sphere and A and a are the constants of integration. 
Whence a particular solution of (6) is 

d ^~^^vt 
xjj^ = At — cos X (r — a + «)• 

tA/t / 

Integrating this with respect to X between the limits co and 0, and then changing A 
into F (a) and integrating the result with respect to a between the same limits, we 
obtain 



'''^~ 2^{vt) dr]o r ®^P' | — 4,,^ 



da. 



4G 
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Performing the differentiation and then integrating by parts, we obtain 



^2 



= -iV 



/.OO 

TT j 
Vt Jo 



F(«) 



+ F^ (a) > exp. i — 



(r — a 4- of 



da 



+ i V ^^ IX^) ^xp. ] — 






00 







We shall presently show that it is possible to determine F(a), so that F(0) = 0^ 
and F(a) e""'*" = when a = oo ; hence the term in square brackets will vanish at both 
Hmits, and we obtain 



v,= 









TT 



exp. 

F(a) 



a' 






T 



+ r M } exp. { - 



(r — a + a)' 



d 



a. 



• ■ (8) 



We must now determine the functions x ^^^ F so as to satisfy the surface conditions 
(2) and (3). 

Equation (2) will be satisfied if 

^(a)-F(a)--.aF'(a) = — . . (9) 



Equation (3) requires that 



Ya^=-i ^^ f ^ X («) exp. ( - £) ^« + 4 V^ f / («) exp. 



a' 



— 7 " 1 da 



a 

2 



\^^ lo ^^ (^) + "^^^ ^"^ ^ i '^^P' ( "" £) ^"- 



Integrating the last term by parts, the preceding equation becomes 



Ya^=i 



Vl \j-^ (") + F («) + aT («) + a^F" (a) } exp. ( - £-) da, . ( 1 o) 



provided, {F (a) + aW (a) ] exp. ( — a^jivt) vanishes at both limits. This requires 
that F(0) = F'(0) =0, and that F(a)€~''' and F'(a)€"'*' should each vanish when 
a = 00 . When this is the case (10) will be satisfied if 



- X (oe) + F (a) + (^' (^) + ^'F^' {a) 



2Va3 



TT 



» ♦ . »; ♦ • \/ 



Whence by (9) 



F"(a) 



3Va 



TT 
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and, therefore. 



F(a) = ^-^' + Ca + D. 



ITT 



The conditions that F (0) = F'(0) = require that C = D = 0; whence 



F (a) = ^^3-, 



iir 



Ya /3aP 



TT 



2 



Also the preceding value of F(a) satisfies the conditions that F(a)€""''', and F'(a)€"~''' 
should each vanish when a = co ; whence all the conditions are satisfied, and we 
finally obtain 



^ = 



2r^{iTvt) 



\ -^ 






3Va sin^ e 



m- + «) 



exp. 



(r — a + a)' 
______ 



da. 



. (12) 



The first integral can be evaluated ; in the second put r — a + a = 2u^/{yt) diXidi 
we obtain 



Ya sin^ 6 I fvt 

xj; z=z — I ^vt + 6a A / — + 



2r 



a 



2 



3V^ sin^ 



V 



TT 



\.,__A-^{'i^U\/vt ~ r 4- af + 2t^\/z/2^ — r + al e"'^" du, (13) 



4. When ^ = the second integral vanishes, whence the initial value of \|r is 



^ = 



Ya? sin^ 6 

2r 



which is the known value of ^ in the case of a frictionless liquid, as ought to be the 



case. 



When t is very large, we may put i = oo in the lower limit of the second integral, 
which then 



ma sin^ e 







1 



Va sin ^6> 

27' 



{3z/« + Qa^ivt) + l(a'^ - r^)}, 



whence 



xh = iYa^ sin^ ^ ( - - - 
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This equation gives the value of xjj after a sufficient time has elapsed for the motion 
to have become steady, and agrees with Professor Stokes's result. 
5. Let Vt be any solution of the partial differential equation 

*(i)«=f (") 



Then, i£ Vq = 0, F(^ — t)v^ dr, where F(t) is any arbitrary function which is inde- 

pendent of r and t, and does not become infinite between the limits, will also be a 
solution of (14) ; for, substituting in (14), the right-hand side becomes 

F{0)Vf + f r{t ^ t)v, dr = F{t)vo + f F(^ - t) P dr 

J Jo ^'^ 



d 



't 



^ y ]/('- '>^ ^^' 



if Vq = 0. 



6. The second expression on the right-hand side of (13) is the value of i/f^ sin^ 6 ; and 
it is easily seen that this expression vanishes when i^ = 0. Hence it follows that the 
expression which is obtained from (13) by changing t into r and V into ¥'{t — r) dr, 
and integrating the result from t to 0, is also a solution of (1). Now, if F(0) = 0, it 
will be found in substituting the above-mentioned expressions in (2) and (3) that F{t) 
is the velocity of the sphere, supposing it to have started from rest ; hence this expres- 
sion gives the current function due to the motion of a sphere which has started from 
rest, and which is moving with variable velocity F(^). 

In order to obtain the equation of motion of the sphere, we must calculate the 
resistance due to the liquid ; but in doing this we may begin by supposing the velocity 
to be uniform^ and perform the above-mentioned operation at a later stage of the 
process. 

If the impressed force is a constant force, such as gravity, which acts in the direction 
of motion of the sphere, and Z is the resistance due to the liquid, it can be shown, as 
in Professor Stokes's paper, that 

(pa cos ^ — p -^ sin^ 9 ) sin 9 d9, 



\ dt / a 



Z = 2iTa 
and that 

dp , ^d^^lr-^ , ^ 



dd ^ dt dr 



where p is the density of the liquid ; also, since 



r 



p cos 9 sin 9d9 = — - 1^ sin^ 9 -~ d9, 

Q J A d(y 
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we obtain 



Z== 



— TTpa 



A. 

dt 



TT 



d^^ 



a ~~ + 2t/>3 ) sin¥ dQ + M'^ 



a 



a? dt\ dr 



+ ^A + Wg, 



a 



where M' is the mass of the liquid displaced. Now, if V were constant, we should 
obtain from (13) 



'^'S'I=~'^('^^+^'' V^+ 



¥' ]> 



and 



whence 



a 



#1 

dr 



+ 2xli,] = - V kvt + 9a /s/"^ + laA 



We must now change t into r, V into F' (t — t) dr, and integrate the result with 
respect to r from t to 0, and we obtain 



Z 



M' d 

a^ dt 



/vt\ 



F'{t — r) nvT + 9a A^ ~jdT + ^Wv + M'^, 



and the equation of motion of the sphere is 



(M + W)v + 



9M' d 



rt 



a^ dt , 







^'(-^)(i- + -v/"K=(M-M> (>.) 



Integrating the definite integral by parts, and remembering that F(0) = 0, the 
result is 







TTT, 



and, diffeientiating with respect to ^, (15) becomes 



(M + IW) V + 



9M 



/ f 



Vv 






x/-^ 



dri = [M--W)g. , . (16) 



Let or be the density of the sphere, and let 



(<y - p)9 

or •+ 2P 
MDCCCLXXXVIIL —A. 



=/; 



9 



H 



hy X = hv^ 



5 ? « « 



(17) 
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then (16) becomes 

v + \v + haAs^/ ^\^~-^-~^dT=j: ...... (18) 

This is the equation of motion of the sphere, from which F (t) or v must be 
determined. 

7. Up to the ])resent time we have supposed the motion to have commenced from 
rest, so that F (0) = 0. Let us now suppose that the sphere was initially projected 
with velocity V. In order to obtain the equation of motion in this case we may 
divide the time, t^ into two intervals, h and t — /^, where ^ is a very small quantity, 
which ultimately vanishes. During the first interval'^ let the sphere move from rest 
under the action of gravity and a very large constant force, which is equal to 
(M + 2M.')X, and then let the large force cease to act. This force must be such as to 
produce a velocity, V, at the end of the interval, h, whence we must have V = X/i, 
i; = Xt ; and, therefore, v = V^//^. Changing f into f-^ X in (18), multiplying by 
e^\ and integrating between the limits t and 0, we obtain 

ve""^ = — Jca \/ - du e^^^ F(u — r) -r + \ Xe'^'du + / e^^'du. . . (19) 

V ttJo Jn ^ ' v^T ' Jo ' "^ Jn ^ ^ 



Now F'(^) is composed of two parts : a large part which depends upon X, and which 
is equal to V/A ; and another part which depends upony, and which we shall continue 
to denote by F'(i^). Hence (19) may be written 



ve^^ = r 



(^^^ - 1) -I- ^{^^ - 1) - Jca \/l\' du |V(t^ - r) € 



AU 



V'^ 

ha A / -- e^''x{u)du, (20) 



where 



/ , [^ Ydr 



o^\/' 
Now X (^0 depends on X, and therefore vanishes when it > h. When u < A, 

X(u) — 2Yu^/h; 



therefore 



e^^X (^) du = -j~ ii^t^^'du = 0, when h = 0. 



ft 




Hence, in the limit when h vanishes, (20) becomes 



* The following procedure, suggested in a Report upon this paper, has been substituted, for the 
remainder of this section as originally written. 
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and the value of the acceleration is 



V—— VXe-^* 



+/e"" - ^ V ^ 



d 



TT dt J 



du 



"ii 





dr 



(22) 



8. It seems almost hopeless to attempt to determine the complete value of F from 
the preceding equations, but, in the case of many liquids, z^ is a small quantity, and 
(22) and (23) may then be solved by the method of successive approximation. For 
a first approximation 

whence 



f 



F^ {t - t) dr 



= / 



^'•* 



s/ii - r) 



(23) 



The integral on the right hand side of (23) cannot be evaluated in finite terms, and 
we shall denote it by ^ (i^). Putting r = ty, we obtain 



where 



Now 



Therefore 



and therefore 



A/y 



dy 






H 



1.3 . . .(2n-- 1) 



n 



T%\ 



Jo 



— Kty 



dy 



1 -6 



\t 



f 

Jo 



yn^-Uy ^y^ ( — ) 



n 



U 



d \«1 — e-^* 



d.\t. 



j, (t) = ^/t 



\t 



Xt 



+ t{- yiin 



\t 



d V I — 6- 



x.t 



d .\t, 



\t 



. (24) 



. . (25) 



When t is very large we may replace (1 — e '^*)/Xi! by {\t) ^, and we shall obtain 



*«=v<{' + ^.w 



\^^(\t - 1)' 

which shows that (f)(t) = where t = co , 

Another expression for (f> (t) may be obtained in the form of a series, for 



*('^='-"L'-^=^^' 



m (2\ty 

H 2 



( - Y (2X0 



n 



» • • 



1.3 , , . (2n-\- 1) 



I • • • 



, (26) 
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by successive integration by parts. The above series is convergent for all values of t, 
and is zero when i^ = cx). 

For a second approximation, (22) gives 



v = ¥{t)=fe-^*-f]ca^l 



d 
dt 



rt 



\u 



(f)(t -- u) du, . . , (27) 



and 



V =: Ye-^ 



+ /(!.-- e-^t) -^fha /^/~ r €~^" (l>{t-u) du, . . (28) 



Let 



x{t) = 



d 



n 



dt . 



e-^"- ^{t — u) du, 



and (27) becomes 



r{t)=fe~-*-fka fixity 



(29) 



Whence to a third approximation 



;= -Y\e-+fe--fka ^^ x(0 +^||/«£^-^''-"'x(« - t) ^^ 



V 



s/'^ 



Let 



m = 



\/{^ - '^) ' 



» • • »• » ».♦• \ tj \j J 



and the last equation becomes 



V 






and 



» » 



• (31) 



/ 



A 



— A< 



/^^ 




P 



IT 







^^^^(^ — 'U)(lt^ + 



fWciH 



IT J 







^''xlj{t -- u)du. (32) 



We must now express all the above integrals in terms of <^(^), From (29) we 
obtain 

d [^ 



X(0 



dt] 







■''^'-''^{u)du 






4> (0 - Xe' 



A^ 



'i p , dr 

du ^' -- 
h v'^ 
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by (24). Changing the order of hitegration, the last integral 






^ . d^c 



At 



^ T 



V- 



-0 



^ X ( ^ 
^r 



=:e-U(OU + 



2\. 



V- 

\ 



— ./r I d\ 



whence 



X (0 = (4 - ^i) i> (0 + Vt. 



• • • > 



(33) 



Substituting this value of x(0 ^^^ i^^)f ^^ obtain 



Now 



also 



and 



whence 



Again. 



xjj{t) 



'(l-Xr)c^(r) ^" 





v/(«-t) + .o 



V 



^-T 



C?1 



^/(^ - t) 



C?1 







Am 



rZii 



o\/{(^-t)('^-'^0} 



cZ'U 



r^; 







e-^"" dr 



u 



y{it-r)(T-u)} 



=z 7T \ € 





— Aw 



TT 



C^M = ^ (1 - e-^'), 



f 






dr 





T6 



Ait 



(^^6 



r?? 



=: d 



U 




TT 

2 

TT 

2\ 










(t + t^) 



■^(l_2e-^*)+.i(l-e-^0}. 



fiv 



t — T 



dr = ijTTt, 



t/;(/5) =: TTi^e-'^^ 



9 * 



« 9 



Jo Jn 







whence (31) and (32) finally become 



• • (34) 



• » ly 



(35) 



• t«9 •»•!> 



(36) 



» • 



(37) 



»«>»«»» 



(38) 
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V =/e~^* — VXe""^* --fka A/ ^ 



V 



(i - X0#) + \^i \ +/&%'^^€-^*(l -~ iX^), (39) 



/ 



{1 



c-^^) + Ve-^* 



-fka aJI 



t + ^A<l>{t) -- ^ I + lfkH^vth-^\ . (40) 



These equations determine to a third approximation the values of the acceleration 
and velocity of the sphere, when it is projected vertically downwards with velocity, 
T, and allowed to descend under the action of gravity. If the sphere is ascending 
the sign of g must be reversed. 

If no forces are in action we must puty= 0, and the preceding equations give the 
values of v and t; to a first approximation only; but, on referring to (21) and (22), it 
will be seen that the values of these quantities to a third approximation may be 
obtained in this case from (39) and (40) by changing/ into — VX and expunging the 
terms/e""^^ and/X"*^ (1 — e""^*). We thus obtain, since X znhv^ 



t; = — Yh 



V€ 



■kt 



rj i I t— 



V a/c^v^ 



{(i - ^ 4>{t) + %/i} - Ya^^pHe-^' (1 - i^t), (4 1 ) 



V 



^ £ f l U l f l W 



Ya¥p^ 



s/tt 



*^+ 0^)^(0 



2X 



X 



^Ya^¥vHh ^\ 



(42) 



9. It appears from the preceding equations that the successive terms are multiplied 
by some power of h as well as of v. If k is not a very large quantity, and the velocity 
of the sphere is not very great, the foregoing equations may be expected to give fairly 
correct results ; but if ^ is a very large quant' ty, it may happen that, notwithstanding 
the smallness of p, kv may be so large that some of the terms neglected may be of 
equal or greater importance than those i*etained. Now, from (17), ^= 9p(2cr + p)~^<^'~^5 
if, therefore, the sphere is considerably denser than the liquid, k will be small provided 
a be not very small ; but if the sphere be considerably less dense than the liquid, k 
will approximate towards the limit 9a~^, and this will be very large if a be small, and kv 
may therefore be large. On the other hand, it should be noticed that when kv or X is 
large the quantities €~^* and ^{t) diminish with great rapidity, and it is therefore by 
no means impossible that the formulae may give a fairly accurate representation of the 
motion even in this case. 

All that we can therefore safely infer is this, that in the case of a sphere ascending 
or descending in a liquid whose kinematic coefficient of viscosity is small compared 
with the radius of the sphere (all quantities being of course referred to the same units), 
the formulae would give approximately correct results, provided the velocity of the 
sphere were not too great. But, in the case of small bodies descending in a highly 
viscous liquid, it is possible that the motion represented by the formulae may he very 
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different from the actual motion ; and if this should turn out to be the fact, the 
solution of (18) applicable to this case must be obtained by some different method. 

Equation (39) shows that after a very long time has elapsed the acceleration 
vanishes, and the motion becomes ultimately steady; in other words, the acceleration 
due to gravity is counterbalanced by the retardation due to the viscosity of the liquid. 
When this state of things has been reached, the terminal velocity of the sphere is 

/ 2^2 /a ^ \ 
X 9v \p /-^ 

This agrees with Professor Stokes's result, who applies it to show that the viscosity 
of the air is sufficient to account for the suspension of the clouds. 

10. We shall now consider the motion of a sphere which is surrounded by an 
infinite liquid, and which is rotating about a fixed diameter. 

We shall begin by supposing that the angular velocity of the sphere is uniform 
and equal to <y, and shall endeavour to obtain an expression for the component 
velocity of the liquid in a plane perpendicular to the axis of rotation, on the supposi- 
tion that no slipping takes place at the surface of the sphere. 

Assuming that the motion of the liquid is stable, it is easily seen that none of the 
quantities can be functions of <j6, where r, 6^ and ^ are polar coordinates referred to the 
centre of the sphere as origin. If, therefore, we neglect squares and products of the 
velocities, the component velocity, v\ of the liquid, perpendicular to any plane con- 
taining the axis of rotation, is determined by the equation 



cU \ dr^ r dr r^ sin d0\ dO) r^ sin^ 6 



>, 



and if in this equation w^e put v^ =^ v sin 9, where t^ is a function of r and t only, the 
equation for v is 

d^v 2 dv 2v 1 dv //toX 

dr'^' r dr r^ v dt' •••*•*•• v / 

The value of the tangential stress per unit of area which opposes the motion of the 

sphere is 

^p / 1 d& dv' v' 

i =: — Z/^ "^'"/j Tj. + 7 

' \r sin C7 dq> dr r 

where E is the radial velocity ; but, since R is not a function of ^, the value of this 
stress depends solely on that of v\ Now Professor Stokes has pointed out that 
unless the motion of the sphere is exceedingly slow, the motion of the liquid will not 
take place in planes perpendicular to the axis of rotation, but the velocity of every 
particle will have a component in the plane containing the particle and this axis. But 
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since this component does not produce any effect on the motion of the sphere, which it 
is our object to determine, we may confine our attention solely to the calculation of v. 
In addition to (43), /y must satisfy the conditions : 

(i.) At the surface of the sphere v =^ ao for all values of t. 

(ii.) When t= 0^ v =^ for all values of r greater than a, the radius of the sphere. 

Let V = Re"^^''^ where R is a function of R alone ; substituting in (43), we obtain 



the solution of which is 



whence 



cm 



•n*OvEa 



2 ^E 2E 



r dr 



1-ii 



+ X^R = 0, 



dr 



- cos \{r — a -\- a) 



V 



A 



d^ 

dr 



r 



.— AVi 



cos X(r — a -{- a) 



Integrating this with respect to X between the limits oo and 0, and then changing 
A into F(a) and integrating the result with respect to a between the same limits, we 

obtain 

(r — a -\- of 



V 



^ V vt dr r Jo ^ ^ ^' 



4:Vf 



da. 



Performing the differentiation and then integrating by parts, we shall obtain 



V 



1 

2r 




.00 



IT 
Vt JO 



F(.) 



r 



+ F' (a) 



> exp. < 



(r — a -{• aY 

4:Vt 



[ da, 



(43a) 



provided F(0) = and F(a)£ *"' = when a = co . 
The surface condition (i.) will be satisfied if 



F(a) + a¥\a) =^ ~ 



2d'^(jd 



IT 



whence 



F(a) 



2a^o} 



7T 



(1 



— a/cA 



the constant of integration being determined so that F(0) = ; this value of F(a) 
also satisfies the condition that F(a)€"~''^ =: when a = cxd . We therefore obtain 



V 



a?(o sin Q 

r\/{iTvt) 



.CO 



a 



+ 1 



a' 

r 



■ — a/a- 



exp. 



r {r — a -j- oLf 



4:Vt 



> da. 



(44) 
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Putting r — a -jr a = 2u^[vt) this becomes 



- = --77^.L_4-+ (l --) exp. (- ^-^^-^ )\.-du. (45) 

If r > a it follows that v = when i^ = 0. When r = a and ^ = the lower 
limit of the definite integral (45) becomes indeterminate ; but since, in this case, we 
are to have v = ao sin 6, it follows that if we put k = r -- a the quantities k and t 
must vanish in such a manner that when k = and ^ = 0, hj2.s/{vt) = 0. 

When t = CO we obtain 

, a^co sin . . 

V = . , . . , . . . . (46) 

This equation gives the value of v after a sufficient time has elapsed for the motion 
to have become steady, and agrees with Professor Stokes's result. 

11. Since the tangential stress per unit of area which opposes the motion of the 
sphere is 

I = — vpa -— - 

dr \r , a 

the opposing couple is 

' Udr\Tla 

dT Y'/a J a 



If, therefore, the sphere be acted upon by a couple, N', its equation of motion will 
be 

iV^'^o, + G = N', 

or 

<Taa) d /v\ ^ , : 

op dr \r/a ' 

where 

When the motion of the sphere commences from rest the value of Ti or v^ cosec Q will 
be obtained from (45) by changing t into r, o> into Y {t —r)di:, and integrating the 
result with respect to r from t to 0, where F {t) is the variable angular velocity of the 
sphere. 
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NoWj 

cl /v\ 1 dv V 



dr \rJa a dr a- 



Hence, if co were uniform we should have 



'dv\ ^ . 2« ''"' 






ao) 



Putting u + ^y {vt)a = /3, the definite integral 



"00 



/•OO 



€ 



.r^/a^ 



V^T v/(x^O , W 






a ' Sa^ 



'sJ'iT \/vt vi s/lT 



« • • 



if ^^ be small ; whence 



^dv\ 2m ( , ^ vt^iA aw 



Changing t into r, and oi into F' (i^ — r) dt, (47) becomes 



aao) 



5p 






Putting 



(48) becomes 



— ^ = k, kv = X. 

era" 



CO + Xo) + "-'— - E'' (^ — t) ( a/t — ;:- a / ^7^ I <^'2" 

a v/tt Jo ^ ^ \^ 2a V / 



+ iX-a/^/^f^F(^-T)^^ = ltoN. (49) 



Now we have supposed the motion to have commenced from rest under the action 
of the couple N' ; but if the sphere had initially been set in rotation with angular 
velocity 11^ and then left to itself, it can be shown in the same manner as in § 7 that 
the equation of motion would be 
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where F (0) = fl. Putting 9 (t) for the last two terms, and integrating, we obtain 



Ce) = Oe""^* 



-f 



K{t-~u) 



6{u)du, ........ (51) 



CO = 



dtt J n 



« « 



• (52) 



For a first approximation we have 



ce) = Oe-^^ c5 = — XOe"^^' = F' (t). 



Whence, if ^, x^ ^^d t/i have the same meanings as in §8, a second approximation 
gives 



• (53) 



0) = Qe"^^ + "^ 



2y/' 



TT 



[e-^<*-"'^(M)cZw (54) 

Jo 



And a third approximation gives 



di = — ^fO€~^* + 






x(0 + 



^^z^* ^ 



rt 



2a^ir dt . 



du 







u 







.—kit-'t) 



y^TC^l 



W'c^v^ d [^ 



4:7f dt 



u\ ^""^^'"^''^^W^^^ • • (55) 






2^7r 



Now we have shown in S8 that 



2a\/nT 

""" 47r Jo 



du e 

Jo 



-A(i-T) 



^/rdr 



.—■K{t—u) 



\}j (u) du. 



Also 



£ €-"^<^--> (^ (^) dti = ^ (0 6 + |~) 



A. 



du €""^^^""^^a/tC?T 
Jo Jo 



C^l 



■X«-t) 



•'r 



^/rdu 



-u 



n 







At 



(Vt-t*)cZt 



-k[m 



^"+^<jf>(oi- 



■ (56) 



And the value of the last integral in (56) is given by (38) ; whence 
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OJ = 0€""^^+ r^ 

^ 2v/7r 



m{t+^^-'^' 



2X X 



which determines the value of the angular velocity as far as v^. 

12. Helmholtz and Piotrowski disco \^ered from their experiments that in the case 
of many liquids slipping takes place at the surface of the solid ; when this happens, 
the surface condition is 

fi{v-aoj) = ppajJ^yj^^ . (58) 

where /8 is the coefficient of sHding friction. Putting k = vp^"^, we obtain from 
(43a) 



dr\T]a ^ 'V vt]^ [ a'^ ' d^ J ^^^-^* \^ 4tvt)^^ 

-iV'j:{^+F'(«)}i«p-(-£)* 

provided F (0) = F' (0) =0, and F [a) e""""' and F' (a) e""""' are zero when a = go 
Equation (58) will be satisfied if 

\a k Vcr ka irlc 



the solution of which is 



F =: — 4- Ae^'* + Be'^^ 



where p and g are the roots of the equation 



3 . 1 



^' + (^ + .-j" + b + .-«) = o- • ^'') 

The roots of (59) will be real if a >^, that is, if a>pp/fi. Now, if there is no 
slipping, /3 will be infinite, and therefore, when there is comparatively little slipping, 
/3 will be large, and this relation will be satisfied unless a is small or v is large ; on 
the other hand, if there were no friction between the surface of the sphere and the 
liquid, /3 would be zero, but it seems improbable that any liquid exists which possesses 
the property of viscosity with regard to the internal motion of its particles, and which 
at the same time is incapable of exerting any action in the nature of friction against 
any surfaces with which it is in contact. If therefore /3 were zero, v would probably 
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also be zero, and the liquid would be frictioiiless. We shall therefore assume that the 
roots of (59) are real. 

The constants A a,nd B must be determined from the condition F (0) = W (0) = 0, 
whence 

F (a) = - J^^- 1 1 + 9^':-^^'] 



F (a) = 



(Sk + a) TT 

, : 2aPoy 
irJc (p — 5') 



UP'' 



€^^^'), 



also this value of F satisfies the conditions that F (a) e"'"', and F' (a) e""*' should 
vanish when a = oo : whence the value of v is 



.00 



V 



a^Q) sin 



r^iirvt) }q 



a 



3 



r (8k + a) 



1 -j~ 



ge 



pa 



pe 



<ia- 



p>-q 



MM*B"» 



e — 6 



^^(P-2')_ 



exp. \ — 



4:Vt 



di 



> aa 



• (60) 



13. We shall lastly consider the motion of liquid contained within a sphere, which 
is rotating about a fixed diameter, when there is no slipping, and when the angular 
velocity is uniform. 

In this case v must satisfy the diflPerential equation (43), and also the condition (i.) 
of § 10 ; but (ii.) becomes '^ = when ^ = for all values of r < a : also we have a 
third condition, viz., that the velocity must be finite at the centre of the sphere. 

A particular solution of (43), subject to the condition of finiteness at the origin, is 



V 



= iA/>/ 



TT d 1 

vt dr r 



exp. 



(r — a)' 

4:Pt 



exp. 



4:Pt 



whence if p and q are any quantities which are independent of r and t, a solution of 
(43) is 



exp. < •— 



(r — oi)^ 
Apt 



- — exp. < 



(r 4- a)^' 



4:Pt 



da 



/it d 1 p^ , , 
•^ \ pt dr T }g^ ^ ^ 

d 1 r r^ 

= ,- - d\\ ¥ (a) €- ^"'' {cos X {r — a) — cos X {r + a)} da. 

If we put p = a, q = 0, F (a) = a, the double integral when ^ = is equal to r by 
Fourier's theorem, for all values of r between a and 0, If we put p = oo^ q = a^ 
the integral when ^ = is zero for all values of r which do not lie betvreen a and oo . 
The solution of the problem is therefore contained in the formula 
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exp. 
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exp. ^ 
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dot, (61) 
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where A is a constant, which, together with the function F (a)^ must be 
so as to satisfy the conditions of the problem. 

14. Though I am convinced that a solution of the problem exists in t 
definite integral, I have not succeeded in obtaining it ; and therefoj 
solution of a different character. 

Let S(r) denote the spherical function d[r~'^ mir)jdr ] then a sokv 
subject to the condition of finiteness at the origin, is 



V ^ tK^-""''^ (Xr) + o>a, 



» 4 



when r ■==-ay v =o)a for all values of t, whence 

S (Xa) — 0, 

and the different values of X are the roots of (63). 
Initially v =^ 0, whence 

(oa = — SA;^S (Xr). 



Let X and [x be different roots of (63), and let T = S (/xr), then, since S 
the equation 



dr^ T dr r^ ' 



we obtain, 



'« 







dr dr 



a 



= 0. 







and since by (63), S and T both vanish where r = a, we obtain 



[^STrMr = 0, . . . 



. » 



provided X and /x are different. To find the value of the integral whe 
jLt = X + <^X ; then from (64) 



2Xc?X 
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#S 



d^d^' 



dr dX dr d\ 



dX= 0, 







or, 



f^Q^r'dr = ia'S'^\a), 



« t s » > 



where the accents denote differentiation with respect to \a \ whence 



i A.a^S's (Xa) = 



wa 



'^ d sin \r 







dr T 



d7\ 



0) . 

- sm Aa. 
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Therefore 

. ^ 2m (sin Xa — Xa) 

and 

2m . e-^^'f'* (sill Xa—Xa) S (Xr) ,^^. 

^=a,a-^S ^^^^^ , ..... (60 

whence the Telocity of the liquid, which is equal to v sin 6, can be found. 

When the angular Telocity is Tariable, the value of the retarding couple, and the 
equation of motion of the sphere, can be obtained by a process analogous to that 
employed in I n. 

[March 10th, 1888. — Since this paper was read, a paper has been published in the 
* Quarterly Journal of Mathematics,'"^ by Mr. Whitehead, in which he attempts to 
deTelope a method of obtaining approximate solutions of problems relating to the 
motion of a tIscous liquid, when the terms iuTolving the squares and products of the 
Telocities are retained ; and he applies his method (see p. 90) to obtain expressions 
for the components in the plane passing through the axis of rotation, of the Telocity 
of a tIscous liquid, which surrounds a sphere which is rotating about a fixed diameter, 
when the motion has become steady. It will be observed, howeTer, that the expressions 
for these components contain the coefficient of viscosity as a factor in the denominator, 
and therefore become infinite when the liquid is frictionless. It would therefore 
appear that the method of approximation adopted is inapplicable to the problem 
considered.] 
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ERRATUM. 



' Phil. Trans./ vol. 179 (1888), A. 
Page 63, for equation (67) as printed read 



r. ^ e ^^"^ sin \a S (\r) 



